Perturbative zero-point energy for a cylinder of elliptical section 
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We examine the Casimir effect for a perfectly conducting cylinder of elliptical section, taking as 
reference the known case of circular section. The zero-point energy of this system is evaluated by 
the mode summation method, using the ellipticity as a perturbation parameter. Mathieu function 
techniques are applied. 
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I. INTRODUCTION 

Vacuum fluctuations of quantum fields caused by the 
presence of boundaries produce changes in the zero-point 
energy, which give rise to the Casimir effect. Interest 
in quantum vacuum manifestations, including this phe- 
nomenon, has been propelled by new theoretical and ex- 
perimental advances 0, 13 • Hypothetical prospects of 
technological applications make even more desirable the 
knowledge of fundamental aspects of the theory, such as 
the value of the vacuum energy or its dependence on any 
of the problem conditions (even the sign of the effect is 
hard to predict, although for interactions between dielec- 
tric bodies some progress has been made in Ref. 3). 

Particularly striking is the modification of the zero- 
point energy caused by a change in boundary shape, 
even at an infinitesimal level. As Casimir energies prove 
to be very sensitive to purely geometrical modifications, 
the subject deserves further consideration. This ques- 
tion has already been addressed for a boundary depart- 
ing from spherical, which has implications for QCD flux 
tube models Q- 

On the other hand, cylindrical boundaries of circu- 
lar section have been object of attention under a vari- 
ety of settings: the perfectly conducting case d H @> 
dielectric media with or without light-velocity conser- 
vation IH El El El El El EI, disper- 

sion [lM 1191 120 | . semitransparent boundary |2lJ, varying 
Robin conditions |22J, finite temperature 1 |24| . coaxial 
surfaces [25j, cosmic strings [2(|, etc. However, devia- 
tions from circular section constitute a largely uncharted 
land. In the present work, we make an incursion into this 
territory by considering a perfectly conducting and in- 
finitely long cylindrical surface of elliptical section, which 
slightly deviates from circular shape 2 . 

In Sec. 2 the solutions to the Maxwell equations in 



1 For a general discussion abo ut t he classical limit in temperature- 
dependent systems s ee Ref. |2.'I . 

2 Yet, there are works 27, 2^| where the elliptical coordinates from 
scattering theory play a role in the study of General Relativity 
(averaged) null energy conditions for a Casimir system. 



terms of Mathieu functions are considered, and the 
boundary conditions which determine the eigenmode set 
are established. The summation of these eigenmodes is 
studied in Sec. 3, thereby obtaining the zero-point energy 
of the electromagnetic field as a perturbative expansion 
in the ellipticity. Sec. 4 offers, as an alternative view, 
the spectrum modification on the basis of a conformal 
transformation relating elliptical and circular cases. The 
conclusions follow in Sec. 5. Relevant properties and re- 
sults about Mathieu functions have been included in an 
appendix. Natural units (ft = c = 1) are used through- 
out. 



II. SOLUTIONS TO MAXWELL'S EQUATIONS 
WITH PERFECT CONDUCTOR BOUNDARY 
CONDITIONS 

We will consider the classic problem: find the zero- 
point energy of an electromagnetic field due to a perfectly 
conducting, neutral surface S. In the absence of any 
charge or current density and with time dependence given 
by exp (— iuit), Maxwell's equations in Heaviside-Lorentz 
units are 



V • E (x) = 0, 

V • B (x) = 0, 

V x E (x) = ilu B (x) , 

V x B (x) = -iwE(x) , 



(1) 
(2) 
(3) 
(4) 



where E (x) is the spatial part of the electric field and 
B (x) is the spatial part of the magnetic field. Both vec- 
tor fields satisfy the vector Helmholtz equation, that is, 



(V 2 + W 2 )A(x) 



0. 



(5) 



where A (x) can be E (x) or B (x) . The vector Laplacian 
is 

V 2 A(x) = VV • A(x) - V x V x A(x) . (6) 
The boundary conditions are 



n x E (x) 
n • B (x) 



0. 

= 0, 



(7) 

(8) 



2 



where n is the outward-pointing, unit normal vector to 
the surface S. In particular, we take S to be an infinitely 
long cylinder, along the z-axis, with an elliptical section 
given (in Cartesian coordinates) by the equation 



2 2 

a 2 6 2 ' 



(9) 



where a is the semimajor axis and b is the semiminor axis. 
We introduce elliptic cylindrical coordinates, (£,77, z), a 
left-handed coordinate system related to Cartesian coor- 
dinates by 



x = f cosh (£) cos (77) , 
V = f sinh (£) sin (77) , 



(10) 
(11) 



where 0<£<oo,0<77<27r and / is the focal length 
given by 



/ = y/a 2 - b 2 . 
The ellipticity is defined to be 

/ 

e = — . 

a 



(12) 



(13) 



To avoid confusion with the exponential function, its full 
form, exp, will be used. In the elliptic cylindrical coordi- 
nate system the ellipse © adopts the form 



Co - cosh" 1 (1/e) 



The boundary conditions become 





-Co) 


= 0, 




-Co) 


= 0, 


B 6 (t 


= £0) 


= 0. 



(14) 



(15) 
(16) 
(17) 



While the vector Helmholtz equation cannot be solved by 
separation of variables in all the orthogonal coordinate 
systems the scalar Helmholtz equation can, the elliptic 
cylindrical coordinate system is one of the cases where it 
can pflj . Furthermore, since the media inside and outside 
S are the same, and thus the speed of light is equal on 
both sides; we may split the solutions into transverse 
electric (TE) and transverse magnetic (TM) modes 0, 
IT4L IT5L Each vector component can be written in 

terms of Mathieu and modified Mathieu functions. 



Interior field 



TM modes (B z = 0) 



It is enough to start considering the solution for one of the nonvanishing field components, which we choose to be 
E z . The general form will be js3 : 



E z (£, 77, z) 



00 » 

Yl / 2~ C » (fc) Mc " } & ?) CG " (r '' 9) 6XP ^ ikz) ' 

Yl J ^ S n ( fc ) Ms n ] (€. 9) Se « (V, q) GXP (tfcz) , 



(18) 



where (Mc^ 1 ** and Ms^ 1 **) ce„ and se„ are the even and 
odd (modified) Mathieu functions of the first kind, re- 
spectively. Notice that there is no zeroth order odd Math- 
ieu function. The coefficients C n (k) and S n (k) are con- 
stants. Equation (|16|l implies that 



Ms« fo,^™) 



(19) 
(20) 



where the superscripts refer to the interior region and TM 
mode, respectively. The second subscript indexes the pth 
zero. The tilde distinguishes the odd eigenmodes from 
the even ones. The dependence of the eigenmodes on the 
ellipticity has been suppressed. The eigenfrequencies are 
obtained through the relation 



4g 



k 2 . 



(21) 



The other two boundary conditions are automatically 
satisfied because E and involve the same vanishing 
factors as in E z [30j. 



2. TE modes (E z = 0) 

The TE mode is obtained by the principle of duality; 
B z will have the same form as equation l|18|) . Since 



d 



equation l|15(l implies that 
d 



B, 



(22) 



(23) 



3 



Therefore, 



III. REGULARIZED ZERO-POINT ENERGY 



Me«'(C ,^ E ) 



MsW'(£ 0) £ 



n,p ) 



(24) 
(25) 



where the prime means differentiation with respect to 
£. The other boundary conditions are also satisfied; it 
follows that vanishes on the boundary because oc 
by definition of a TE mode. 



E v and E 



B. Exterior field 

Although this part of the solution is not discussed 
in |30l ] , it can be worked out in terms of functions akin to 
Hankel functions. The reason for doing so is the physi- 
cal demand that, at large distances, the field components 
should behave like cylindrical waves. This is easily met 
by choosing the linear combinations 

Mci 3 > (£, q) = Mc« (£, q) + i Mcf (£, q) , (26) 
M(£> (f , g) = Mel 1 ' (£, g) - i Mcl 2 ) (e, g) , (27) 

with corresponding relations for the odd modified Math- 
ieu functions. 

There are many parallels between modified Mathieu 
and Bessel functions. As with Bessel functions of the 
second kind, the modified Mathieu functions of the sec- 
ond kind are not regular in the interior region. This is 
why they were not included in the solution for the field 
inside the cylinder, but shall now be used for the field 
outside. As we want the waves to be outgoing, the mod- 
ified Mathieu functions of the third kind are selected. 



1. TM modes (B z = 0) 

The z component of the electric field will have the same 
form as equation I jlSj l but with modified Mathieu func- 
tions of the third kind. The boundary conditions imply 
that 

Mc! 3 ) 0, (28) 

Ms^ ,g™)=°> (29) 

where the first superscript refers to the exterior region. 

2. TE modes (E z =0) 
Analogous to the interior TE modes, we find 



Mcl 3 )'fe,g.^ TE )=0, 

M S i 3 )'(e ,Cp TB ) = °- 



(30) 
(31) 



In natural units, the zero-point energy per lateral unit 
length amounts to half the sum (including integration 
for k) of all the eigenfrequencies given by equation 121(1 . 
Since such a quantity is divergent, some form of regular- 
ization is called for; we use zeta-function regularization. 

The even zeta-functions are 



°° pOO 77 



where A € {I, II} and B e {TE, TM}. Let 
4g 



P 



then equation (|21|l becomes 



7 



ui = \/"f 2 + k 2 . 
Integrating over k leaves 



^=27 B 12' 2 



(32) 



(33) 



(34) 



(35) 



P =i 



where B is the beta function. Following similar analysis, 
the odd zeta-functions are 



1 s- 1 
2' 2 



(36) 



There is no odd zeta-function for n = 0. The regularized 
zero-point energy per lateral unit length as a function of 
ellipticity is 



£ c {e)= lim \ Y, 

A.B 



(37) 



where the limit means analytical continuation to s = — 1 . 
As we are summing over the interior and exterior re- 
gions we expect the analytical continuation to be finite 
(see an d references therein). For this reason, nei- 

ther arbitrary scales nor additional prescriptions have 
been introduced. 

In terms of the dimcnsionless variable z — aj, the even 
zeta-functions with .4 = 1 and B = TM are 



p=l 



2tt 



The summation over p can be written as a contour 
inte gral with the help of the argument principle (see 
e.g. lH Hill) as follows 



4 



gK.r]'-^/ an — 



where, using equations (|13|l and l|33|) . 



9 = 



(39) 



(40) 



The integration circuit is the boundary of a region Q of the complex z-plane which contains all the wanted zeros and 
avoids the origin. Other A and B and the odd counterparts follow in a similar fashion. 



A. Ellipticity expansion 



The modified Mathieu functions appearing in the contour integrals can be expanded for small ellipticity in a formal 
series (see Appendix II B|) . Equation l|39|) becomes 



^ [ Z n,P 

■p=i 



TM 



1-s S — 1 



27r * Jdl! 

The corresponding odd sum is 

I l-s S — 1 



dz . 



ln(J„ (*)) 



4J„ (z) 



2?ri 



z 



on 



(41) 



(42) 



which is identical (up to O (e 4 ) ) except for the sign in front of the Kronecker delta. If the even and odd sums are 
added, then 



EK M ]"+EK M ] 
p— i p— i 



= d (n) / dz z 

2tt« Jon 



where to take into account that there is no odd sum for n = 

d{n) 



1, n = 0, 

2, n > 1. 



Integrating termwise (formally), the second term is proportional to the first since 

,zJ' n (z) 



dzz- 

90 J n ( z ) 



(s-1) I dzz- s \n(J n (z)) . 
an 



Thus, equation (|43|) simplifies to 



„_i „_i Z7rz Jao 



1 1— s _|_ r. i i m - 
p— 1 p— 1 

Using equation (|38|) and its odd counterpart, 



dzz s ln(J n (z)) 



s-1 



£Cn™(s) + ££'™( s ) 



n=0 



2tt V 2 ' 2 

n=0 



2th jan 
I — 



dz z s In (J n (z)) 



s- 1 



(43) 



(44) 



(45) 



(46) 



(47) 



The prefactor is exactly what one gets in the circular cylindrical situation. The other combinations of A and 



5 



B all have the same form as equation l|47|l . but with 
their corresponding circular cylindrical prefactors. Thus, 
equation l|3*T|) gives 



£ c (e) = S C (0) 



1 + ~ e z 



0i 



(48) 



where £c (0) is the regularized zero-point energy of per 
lateral unit length of a circular cylinder of radius a. 
While it is doubtful that equation (|48|) converges for 
< e < 1, at worst it is an asymptotic series as e — > 0. 
In either case we may write the next term as O (e 4 ) . 
The numerical value of £c (0) is 



fed.,—*™*" 

which, together with equation (j48(l gives 



£c (e) 



0.01356 



1 + 



O 



(49) 



(50) 



where a is the semimajor axis. Let R = (a + b) /2, then 



2R 



i + VT - ^' 



and 



£c(e) 



0.01356 



1 + (e 4 )] 



(51) 



(52) 



That is, the zero-point energy per lateral unit length in 
terms of R is the same as that for a circular cylinder 
with radius R, up to quartic corrections in the elliptic- 
ity 3 . Equation (|52[) can be related to the existence of 
a conformal mapping in the complex plane which trans- 
forms the ellipse into a circle. This subject is discussed 
in the next section. 



IV. CONFORMAL TRANSFORMATION 



A. Formulation 



Some of the ideas in Refs. 36, 37, 38] suggest the use of 
an adequate conformal map. We shall employ a transfor- 
mation taking the interior of the ellipse 10 to the interior 
of the circle w = Rexp (if), where 



R 



The required map w 



2 

! is |3£ 



' (z) = fiVfcsn 



2K(fc) . 



Va 2 - b 2 



(53) 



. (54) 



3 This is not unique. For instance, R = (a 2 + b 2 ) /2 would cause 
the same effect as R = (a + b) /2, since they only differ at O (e 4 ) . 



where sn is a Jacobi elliptic function, K is the complete 
elliptic integral of the first kind and k depends on the 
semiaxes through theta functions as follows: 



where 



^2 (0,g) 

tf 3 (o,s) 



a — b 



a + b 



(55) 



(56) 



The k and q variables used here should not be confused 
with those in other sections. The theta functions are 
given by 

oo 

tf 2 (*, q) = 2q 1 ' i 1 n{n+1) cos i( 2n + 1) z ) » (57) 
n=Q 
oo 

q n cos (2nz) . (58) 



In terms of e and R, 



(1 = 



= (-)' 
\2RJ 



(59) 



Using the definitions of the theta functions, for small 
ellipticity 4 , 



k = ^ + 0(e*) 



With the help of 



K(fc) = ~ + ^fc 2 + 0(fc 4 ) 
sn (u, k) = sin (u) 



(60) 



(61) 



[u - sin (u) cos (u)} cos (u) fc2 | Q ^ ^ 



the small ellipticity expansion of equation (|54|l is 



w (z) = z - 



Therefore, 



AR 2 



3z 2 



+ 0(e 4 ) 



(63) 



(64) 



where the prime denotes differentiation with respect to 
z. For our purposes, we also need the leading term of 



4 In the opposite case one may wish to consider, e.g., a — > oo and 
fixed b. Using theta function properties, for this limit, w(z) = 
fitanh (ttz/ (46)), which is the transformation taking the region 
between the parallel lines y = ±ib to the interior of the circle. Of 
course, these lines can be viewed as the 2-D projection of parallel 
plates. 
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the inverse of w, given by z (w) — w + O (e 2 ) . With w = 
r exp (i(f), z = r exp (iip) + O (e 2 ) . Using equation l|64[) . 

K| 2 = l-^^ e 2 + 0(e 4 ). (65) 

Under any conformal mapping, the 2-D Laplacian oper- 
ator transforms as 

V* = K| a V». (66) 
We start from the 2-D Hemlholtz equation in the w-plane 



(67) 



where the chosen eigenvalue symbol comes from equa- 
tion l|33|) . The boundary conditions are set on the image 
of the ellipse, which is w — i?exp (i(p), < ip < 2tt. 
Since the boundary is circular, the eigenfunctions are of 
the form 



fn ( r ) ex P ( ini P) ■ 



(68) 



where /„ (r) is the suitably normalized radial function 
which satisfies the boundary conditions. 



B. Rayleigh-Schrodinger expansion 

Next, Dirac notation shall be temporarily adopted 
in order to take advantage of the Rayleigh-Schrodinger 
method. Supposing that the eigenvalues A n and eigen- 
states \U n ) for a Strum-Liouville problem 



L\U n )+X n \U n )=0 



(69) 



with given boundary conditions are known, one wonders 
which are the new eigenvalues A„ and eigenstates \hi n ) 
for the 'perturbed problem' 



L \U n ) + \ n \U n ) - EpcrtVrtl^n) = 



(70) 



under the same type of boundary condition, where £ pcrt 
is some small parameter. Including just modifications to 



first order in £ port , we make the ansatze 



\Un) = Pn) 



K + ■ 



A, 



£pcrtl V n>+0(£pcrt) 
pcrtMn @ ( £ pert) ' 



(71) 
(72) 



and replace them into equation (|70|) . At O (£p 0rt ) equa- 
tion (|69|l is recovered, while the O (£p Crt ) contribution 
yields 



(L + X n I)\V n ) - (r pert -/i„) K), 



(73) 



where / is the identity operator. After applying (U n \ 
on both sides, taking into account the adjoint of equa- 
tion l)69|l (with = L understood), we obtain 



Mr, = (^nkpertl^n)i 



(74) 
(75) 



where dg denotes the integration measure making the 
{U n } orthonormal. 

For our studied case, in the notation of equation l(6"9"|l . 



L = V 2 , 

„2 



(76) 
(77) 



and the integration in equation (|75|l will be on the type 



2-rr p 

dg = I dip dr r, 
o JK 



(78) 



where TZ is the radial range. Now, application of the 
inverse of transformation l|54|) to equation (|69|l , and using 
equation 1)66(1 . give 



V 2 z U n +f n \w'\ 2 U n = 0. 



(79) 



Since the small parameter in the pertubed equation l|70() 
is identified as 



£ pert e i 



(80) 



^2 |„,,/|2 



comparison of equations l|70() and l|79() leads to 7 2 
\i ~ £ port r pcrt; which, taking into account equa- 
tions l(o4 [) . ifTo Jl . (fTT ft and fSD ft . yields 



_ 2 3r 2 cos (2(p) 

'pert — In 2R? 



(81) 



With this r port , we calculate the O (£ 2 crt ) contribution 
to the nth eigenvalue using equation l|75l) 



fj, n <x dip cos (2<p) 
Jo 

= 0. 



(82) 
(83) 



Hence, we conclude that, up to O (e 4 ) , the interior eigen- 
values will not change, and neither will the zero-point 
energy. 



V. CONCLUSIONS 

The main result in this work is (|52p. which shows that 
the Casimir energy per lateral unit length for an elliptical 
cylinder has the same value, up to quartic corrections in 
ellipticity, as for a circular cylinder with radius equal to 
the mean of the two semiaxes. A quadratic correction 
appears if the same energy is expressed in terms of one 
of the semiaxes. 

This can be envisaged from a conformal transforma- 
tion (|54l) which maps the ellipse onto the circle in ques- 
tion and, perturbatively speaking, yields no quadratic 
contribution. Such an infinitesimal symmetry preserves 
the Casimir energy if deformations from circular to ellip- 
tical sections (or vice versa) do not go beyond O (e 2 ) . 



7 



Furthermore, there is a significant connection with Kvitsinsky's work [38J. The 2-D elliptical zeta-function is 



C 2 - D (a, e )=a^ 



,p—i p—i 



(84) 



where a = s — 1. Following a similar analysis as in Sec. IIIII 



O e 4 



where £ 2 ~ D (a, 0) is the 2-D zeta-function for a circular boundary of radius a. If a = 1 + e 2 /2 and 6=1, then 



C 2 " D {o,e) = e~ U (M) 



2-D 



(7 o 

1+ 4 e 



e 4 



(85) 



(86) 



where the 2-D circular zeta-function now has unit radius. Comparison of the different notations leads to e 2 2a, 
a <^ 2p, £ 2_D (f, e) <^> C(Pi ellipse) and £ 2 ~ D (u, 0) <^ ((p;D), where the objects on the right are those in Ref. [38| . 
Then, our result coincides with the unnumbered formula below equation (23) of the referred paper, originally derived 
for p — 2,3,... and just field modes like in our 'I, TM' subset, suggesting that the expression in question could be 
valid beyond its initial settings. 

From the viewpoint of Mathieu functions, one may argue that the obtained relation has been established by virtue of 
a variable change (|53(l in formulas Ij37(l . Unfortunately, the next order in ellipticity is significantly more complicated. 



APPENDIX A: MATHIEU FUNCTIONS 



1. Mathieu functions 



Mathieu's differential equation in canonical form is 

-_d 2 

dr/ 



, , a(q) — 2gcos (2rj) 



H (77) = 0. 



(Al) 



Requiring the solutions to be periodic leads to characteristic values for a (q) . Periodic solutions that are even with 
respect to rj have characteristic values a n (q), whereas those that are odd have characteristic values b n (q). If q = 0, 
then a n (0) = b n (0) = n 2 and the even and odd solutions are the cosine and sine functions. The even and odd 
(periodic) Mathieu functions are, respectively 



ce n (r],q)= ^2 c n.k (<?) cos i( n + 2fc) v] > n = 0,1,2,..., 

k— — 00 
00 

se„(?7,g)= ^2 £ n,k (?) sin [( n + 2fc ) V] > n = 1,2,3, ... . 



(A2) 
(A3) 



k— — 00 



They are normalized such that 



27T C-2TX 

dri[ce n (ri,q)\ 2 = I dr) [se„ (rj, q)} 2 = tt. (4) 
Jo 



The normalized coefficients can be expanded in formal 
power series. For n > 2, the even coefficients are 



„ (g) = l + O(g 2 ), 



c n,±l (l) = T 



1 



4 (nil) 



q + O (q 2 



(5) 
(6) 



For n = 1, we have 



C lfi (q) = l + {q 2 ) , 
Cl,+l(<7)= - \q + 0{q 2 ) 



and for n = 



c Q .o(q) = ^+0{q 2 ) 1 
c 0,+l (q)= -^q + O (q 2 ) 



(7) 
(8) 

(9) 
(10) 
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All other coefficients are O (g 2 ) or higher. The odd co- 
efficients are not the same as the even; however, they do 
agree up to O (g 2 ) with one exception: 

c 2i _! (?) = O (g 2 ) . (11) 



I 

B. Modified Mathieu functions 



The modified Mathieu differential equation is 

d 2 



, a{q) + 2qcosh(2£) 



3(0=0, 



Let £ = 2yfq cosh (£), then 



(C 2 -4q)^ + C^- a (g)-2 g + C 2 



a (C) = o. 



(12) 



(13) 



When q = 0, a (0) = a n (0) = b n (0) = n 2 and equation (|13|1 reduces to the Bessel differential equation. The even and 
odd modified Mathieu functions of the first kind are, respectively 



M<#> (£,«)= d «, fe (9)Jn+2fe(2^cosh(0), n = 0,l,2,... ) 

A;— — oc> 

oo 

Ms^ (£, g) = tanh (£) £ d n<k (q) J n+2k (2^9 cosh (0) , n=l,2,3, 



k= — oo 



(14) 
(15) 



The coefficients are related to those in Section TA II 

d n,k (<?) = P 0?) (-l) fe c n}t (q) , (16) 
dn,k (<z) = P (<?) (-1)" (n + 2fc) c nM (g) , (17) 



Mc« (£,g) C ~°°M S W (£,g) s ~~J n k). 
The normalization factors are therefore 



p(?) 



= C n,k(<l)' 



(18) 



(19) 



fe=-c 
oo 



= J2 {n+2k)c ntk {q). 



k— — oo 



The expressions in equation (|14H and l|15|) are absolutely 
convergence for | cosh(£) | > 1 [30j . 

Using the results from Section IA li the coefficients can 
be expanded in formal power series. The even coefficients 



for n > 2 are 



<o (g) = l~ 2 (n 2 -l) q + °^ 



d n,±i (g) = ± 



where p (g) and p (g) are normalization factors. The mod- 
ified Mathieu functions of the first kind are normalized For n = 1 
to have the same asymptotic form as the Bessel functions 



4 (nil) 



q + O (g 2 ) 



di,o(?) = l+g9 + ^(? 2 ) 
d h+1 (q) = lq + 0(q 2 ), 



and for n = 



^o,o (g) = 1 + ^ 9 + ^ (? 2 ) 

<Vu {q) = \q + o (g 2 ) . 



(21) 
(22) 

(23) 
(24) 

(25) 
(26) 



(20) All the other even coefficients are O (g 2 ) or higher. The 
odd coefficients for n > 2 are 



d n,a (g) = 1 + 

d n ,±x (g) = ± 



1 



2 (n 2 - 1) 



g + 0(g 2 ), (27) 



n± 2 



4n (nil) 



g + O (g 2 



(28) 
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For n = 1 All the other odd coefficients are O (g 2 ) or higher. 

d h0 (q) = l + lq + O(q 2 ), (29) 
d h+1 ( q ) = ^q + 0(q 2 ). (30) 

I 



1. Ellipticity expansion 

Consider the even modified Mathieu functions of the first kind. With £ = cosh -1 (f/e) and q = z 2 e 2 /4, equa- 
tion JHJ| is 

oo 

Mc«(C ,g) = E d nM (^ 2 e 2 /4) J n+2k (z) . (31) 

k— — oo 

For n > 2, using equations lj21|) and 122(1 . the formal series in ellipticity is 

Mc« (Co, g) = J„ (*) + (ra ' 1)J " +2(z) ; 6 2 ( J ; 2 ( ! ) i ; (n + 1)J "- 2(z) Z 2 e 2 + O (e 4 ) . (32) 
Using Bessel function recurrence relations, equation l|32l) simplifies to 

Mc« ( ?0 , g) = J„ (z) - e ' + ° ( e4 ) ' ( 33 ) 

Similarly, for n = 0, using equations (|2*5)l and lt2*6l 

Mc« (Co ,g)=J (z)-i4^ e2 + ( e4 )' 



4 

that is, the n — case follows the same rule as for n > 2. Using equations l|23(l and 1)24(1 ■ the n = 1 case is 



(34) 



(35) 



Mc« (C„, g) = J x (z) - J (ji (z) - i Jo (*)) e 2 + O (e 4 ) . 
Collecting all the results, in general, for n > 

Mc« (Co, 9) = J n (z) - \ (V n (z) - Jo (z)) e 2 + (e 4 ) . (36) 
By similar analysis, the ellipticity expansion for the odd modified Mathieu functions of the first kind is 

Ms« (Co, ?) = J„ (*) - | (j; (*) + S ~f Jo (*)) e 2 + O (e 4 ) . (37) 
The same procedure is repeated for the modified Mathieu functions of the third kind and for the required derivatives. 
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